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Can Spinor Dipolar Effects be Observed in Bose-Einstein Condensates? 
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Weak dipolar effects in atomic Bose-Einstein condensates (BECs) have recently been predicted 
to develop spin textures. However, observation of these effects requires magnetic field as low as 
~ 10 /iG for spin-1 alkali BECs, so that they are not washed out by the Zeeman effect. We present 
a scheme to observe the magnetic dipole-dipole interaction in alkali BECs under a realistic magnetic 
field of ~ 100 mG. Our scheme enables us to extract genuine dipolar effects and should apply also 
to 52 Cr BECs. 



PACS numbers: 03.75.Mn,67.57.Fg,03.75.Kk 



One of the salient features of gaseous Bose-Einstein 
condensates (BECs) is the magnetic dipole-dipole in- 
teraction, which is long-range, anisotropic, and exerts 
a tensor force. The dipolar interaction is expected to 
yield rich phenomena when combined with spin degrees 
of freedom, such as the Einstein-de Haas effect N, [2| 
and ground-state spin textures and mass currents [3|, [4[ , 
because dipole-induced spin textures can generate mass 
transport via spin-gauge symmetry [5j]. Unfortunately, 
however, the observation of these phenomena requires 
extremely low, if not unrealistic, magnetic field typically 
below 10 jiG for spin-1 alkali BECs, so that they are not 
masked by the Zeeman effect. The only dipolar effect ob- 
served to date is the anisotropically distorted expansion 
of spin-polarized 52 Cr BEC 0. 

In this Letter, we show that a weak magnetic dipole- 
dipole interaction of atomic BECs induces spin textures 
even in a magnetic field as high as 100 mG. At such 
high magnetic field, the total magnetization of the sys- 
tem along the magnetic field is fixed by the Zeeman en- 
ergy, and therefore the predicted dipolar effects, which 
rely on the spin relaxation, are unobservable []], Q, 0, Q| • 
However, we show that the dipolar interaction can flip 
local spins within the Hilbert subspace of a fixed total 
magnetization, allowing spin textures to be generated. 
This spin texture can be observed by spin-sensitive in 
situ imaging techniques [7[ or by Stern- Gerlach separa- 
tion. 

We consider the dipole-dipole interaction between two 
magnetic dipole moments fii = gji^F^ and fi 2 = 
gfi^Frf , where F^ h and represent the hyperfine 

spin operators in the laboratory frame, g is the Lande 
^-factor of the hyperfine spin, and /iB is the Bohr mag- 
neton. The dipolar interaction between atoms located at 
ri and r 2 is given by 
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where r 12 = n - r 2 , e 12 = r 12 /\r 12 \, and c dd = 
Mo(#Mb) 2 /(47t) with fiQ being the magnetic permeability 
of vacuum. For spin-1 alkali atoms, we have g = 1/2 [4|]. 

Suppose that a homogeneous external magnetic field 
B is applied in the z direction. Since the linear Zeeman 
term —fi z B rotates the atomic spin around the z axis at 
the Larmor frequency cjl = g^Bjh, it is convenient to 
describe the system in the rotating frame of reference in 
spin space. The spin operators Fj (j = 1,2) in the ro- 
tating frame are related to those in the laboratory frame 
by F j± = F jx ± iFjy = e ±iuj ^F)f and F jz = F)f. The 

dipolar interaction (pp) can be rewritten in terms of F as 
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where Y 2m (ei 2 ) represents the m-th component of the 
rank-2 spherical harmonics. When the Zeeman energy 
is much larger than the dipolar interaction energy, i.e., 
Huj-l ^> Cdd | ^| 2 / 1^12 1 3 , the spin dynamics due to the dipo- 
lar interaction is much slower than the Larmor preces- 
sion. In a large magnetic field, we may therefore use an 
effective dipolar interaction which is time-averaged over 
the Larmor precession period 27t/o;l, 
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In contrast to Eq. (pp), the time- averaged interaction (j3]) 
separately conserves the projected total spin angular mo- 
mentum and the projected relative orbital angular mo- 
mentum. However, the long-range and anisotropic nature 
of the dipolar interaction is still maintained in Eq. ([3]). 
We will show below that (vdd) generates spin texture 
even in the case of weak magnetic dipole coupling of al- 
kali atoms. 

We study the spin dynamics in a system of N spin- 
1 Bose-Einstein condensed atoms with mass M confined 
in a spin-independent potential Ut rap {r) in a magnetic 
field B. In the rotating frame of reference in spin space, 
the dynamics of the system is governed by the non-local 
Gross-Pitaevskii equation 
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where ^ m {r,i) is the condensate wave function for the 
magnetic sublevel m (which is related to that in the 
laboratory frame by ^ b (r,t) = ^ m (f, t) exp^mcc^)), 
Ho = -h 2 V 2 /(2M) + /7 trap (r), and q = (ji B B) 2 / (4E M ) 
is the quadratic Zeeman energy with E^f being the hy- 
perfine energy splitting. We define the number den- 
sity n = |VVn| 2 , longitudinal magnetization f z = 



| , and transverse magnetization /+ = /* 



fx + ify — V%(^t^o + i/jqi/j-i). The short-range inter- 
actions are characterized by go = Ai:h 2 {ao + 2a2)/(3M) 
and gi = 47rft 2 (a2 — ao)/(3M) with a s (s = 0,2) being 
the 5-wave scattering length for the scattering channel 
with total spin s. We note that the linear Zeeman term 
is canceled in the rotating frame. 

The time- averaged dipolar interaction (-Odd) induces an 
effective dipolar field b{r) given by 



b{r) =C dr 



j ^2 (e) 

_ -f/13 
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where C = V&rcad/iVSg^B), e = (r — r')/\r — r'|, and 
6 ± = 6 X ± zfe^. It can easily be verified from Eqs. (|4|) and 
([5]) that the inhomogeneous Larmor precession around 
b{r) changes the local longitudinal magnetization f z {r) 
as well as the local transverse magnetization /+(r), while 
the total longitudinal magnetization J drf z (r) is con- 
served. 

We perform numerical simulations of the spin dynam- 
ics using the parameters of the Berkeley experiment . 
In Ref. [7[, a BEC of spin-1 8T Rb atoms was prepared in 
the m = 1 state. Then the spin was rotated to a trans- 
verse direction by applying a tt/2 pulse of the rf field and 
the ensuing Larmor precession of the transverse magneti- 
zation was probed by spin-sensitive imaging. The crucial 



observation to be discussed in detail below is that a sig- 
nificant dipolar effect can be observed only if the spins 
are rotated through an angle other than tt/2. 

We consider a spin-1 87 Rb BEC of N = 4x 10 6 atoms in 
a trap with frequencies {uo Xl uo yi uo z ) = 27r(150, 400, 4) Hz 
in a magnetic field of B = 54 mG. We first prepare a spin- 
polarized stationary state in i/ji by solving Eq. (j4j) with 
the imaginary-time propagation method while keeping 
= = 0> an d then rotate the spin state about 
the y axis by a tilt angle 6 at time t = 0. We simulate 
the ensuing spin dynamics by solving Eq. (|4j) in three 
dimensions with the Crank-Nicolson scheme. The dipolar 
term is calculated by the convolution theorem with a fast 
Fourier transform algorithm [8[. 

The Thomas-Fermi (TF) radii of the BEC are 
(R x1 R yi R z ) = (6.7,2.5,252) fim and the peak density 
is no = 5.6 x 10 14 cm -3 , which gives a dipole healing 
length of £dd = ft/V2Mcdd^o = 5.7 /mi. Since R x and 
Ry are comparable to £dd, the spin texture discussed be- 
low is almost one-dimensional in the z direction [4|. We 
characterize the spin texture in terms of the integrated 
magnetization f(z) = J J dxdyf(r). 

Figure 1 shows the dynamics of texture formation for 
= 40°. A schematic representation of spin texture is 
shown in Fig. 1(a). Figures l(b)-(e) show snapshots of 
the spin textures, while Fig. 1(f) show the spin texture 
at t = 637 ms calculated without the dipolar interaction. 
Since no texture is developed in Fig. 1(f), we conclude 
that the spin textures seen in Figs. l(c)-(e) arise from the 
dipolar interaction. The helix observed in Figs. l(c)-(e) 
is symmetric with respect to z = 0, i.e., it is right-handed 
for z > and left-handed for z < 0. A comparison of 
Fig. 1(f) with Fig. 1(b) shows that the quadratic Zeeman 
effect slightly rotates the spins. Figure 1(g) shows the 
direction of the transverse magnetization (j) = arg(/+) as 
a function of z. 

To evaluate the development of the spin texture, we 
define the phase difference A<fi = 0(0 /im) — 0(155 jam) 
between z = jam and z = 155 /im, and plot the time 
development of A0 for several tilt angles in Fig. 2(a). 
Initially, |A0| increases linearly and the spin texture de- 
velops more prominently for smaller 6. Since the second 
term in Eq. ([5|) is parallel to f(r) for uniformly polar- 
ized systems, it does not initially contribute to the spin 
dynamics. Therefore, the time dependence of at t = 
can be derived from Eqs. (|4]) and ([5]) as 
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and we neglect the quadratic Zeeman effect. We plot the 
z dependence of bo{z) in Fig. 2(b). The time dependence 
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FIG. 1: (Color) (a) Schematic representation of spin texture. 
The color of each arrow represents the azimuthal angle <\>{z) — 
arg[/+(z)] (see the color legend), (b)-(e) Snapshots of spin 
texture seen from the y direction at (b) ms, (c) 318 ms, 
(d) 531 ms, and (e) 637 ms after an rf pulse of = 40°. 
(f) Snapshot of spin texture at 637 ms calculated without 
the dipolar interaction. The arrows show the projection of 
f(z) onto the x-z plane. The spin texture is displayed from 
z = —214 /im to z = 214 /im. (g) Time development of <p(z). 
Note that <j>{z) at z ~ becomes flat as time develops (see 
text) . 



of A(j) at t ~ agrees well with Eq. (|6|) as shown in 
the inset of Fig. 2(a), where the black line corresponds 
to A(j)/cos6 = —(gfiB/h)Abot with Abo = ^o(0 /im) — 
6o(155 /im), as indicated in Fig. 2(b). 

It follows from Eq. ([6]) that the dipole-induced spin 
texture does not appear when = 90°, indicating that 
we can probe the inhomogeneity of the external field by 
setting = 90°. Thus we can experimentally distinguish 
the genuine dipolar effect from the one arising from resid- 
ual magnetic-field gradient by changing 0. 

As the spin texture further develops, the effective field 
cannot be approximated with bo(z). In fact, A(f)/cos6 
deviates from the asymptotic line at t ~ h/e<id = 260 ms, 




100 200 300 400 500 600 

t [ms] 




-300 -200 



-100 100 200 300 



(C) 1 




-300 



-200 



-100 
z [^m] 



100 



200 



300 



FIG. 2: (Color) (a) Time development of the phase difference 
Acj) between z — /im and z — 155 /im. The inset shows the 
same quantity scaled with cos 0. The black line corresponds to 
A(j>/ cos 6 — —(gnBAbo/h)t. (b) z dependence of the effective 
dipolar field at t = defined by Eq. (0. (c) Profiles of the 
longitudinal magnetization f z (z) at t = ms (dotted curves) 
and t — 531 ms (solid curves) scaled with the number density 
n(z = 0, t = 0) = // dxdyn(x, y, z = 0, t = 0) . 



where edd = 27TCdd^o/3 is the dipolar energy. After 
t ~ h/ead: the spins tend to align in a high density re- 
gion and the helical structures appear only in low density 
regions, so as to reduce the kinetic energy cost [see the 
profile of (j)(z) at t = 637 ms in Fig. 1(g)]. Moreover, 
to reduce the energy cost of the helical structure, the 
transverse magnetization decreases in the region where 
\d(j)/dz\ is large, and instead the longitudinal magnetiza- 
tion increases in this region. In Fig. 2(c), we plot the 
distribution of the longitudinal magnetization f z (z) at 
t = ms (dotted curves) and t = 531 ms (solid curves) 
for several values of 0. We note that the distribution of 
f z (z) is significantly deformed by the dipolar interaction, 
while J dzf z (z) is conserved. After t ~ 600 ms the spin 
texture develops into a complicated structure. 

We have also simulated the texture formation in a 
pancake-shaped trap. We consider a spin-1 87 Rb BEC of 
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FIG. 3: (Color) Snapshots of spin texture in a pancake-shaped 
trap at t — 749 ms after an rf pulse of — 60° . The magnetic 
field is (a) parallel and (b) perpendicular to the symmetry 
axis. The color shows the direction of the transverse magne- 
tization cj) = arg(/ + ), where the overbar denotes an integrated 
value along the symmetry axis. The white circles show the 
TF radius and black arrows show transverse magnetization 
(f x , f y )/n, where the black arrows in (b) represent (f x , f y )/n 
on the x-z plane of the coordinate space with — > and f indi- 
cating (f x , f y )/n — (1, 0) and (0, 1), respectively. 



N = 5 x 10 6 atoms in a magnetic field of B = 50 mG in 
a pancake-shaped trap with (u±,lu\\) = 27r(4.3, 142) Hz, 
where co± and cjn, respectively, denote the trap frequen- 
cies perpendicular and parallel to the symmetry axis of 
the trap. The TF radii of the BEC are (R±,R\\) = 
(100, 3.0) fim and the peak density is no = 1 x 10 14 cm -3 , 
giving R\\ < £dd = 13 fim < R±. Figure 3 shows snap- 
shots of the spin texture at t = 749 ms after rotation of 
the atomic spin about the y axis by = 60°. The symme- 
try axis of the trap is taken to be parallel [Fig. 3(a)] and 
perpendicular [Fig. 3(b)] to the magnetic field (B \\ z). 
Because of the anisotropy of the dipolar interaction, the 
obtained spin texture strongly depends on the direction 
of the trap anisotropy. When the magnetic field is par- 
allel to the symmetry axis, the spin texture is concentric 
in the x-y plane [Fig. 3(a)]. On the other hand, when 
the magnetic field is applied perpendicular to the symme- 
try axis, the spin texture is highly anisotropic [Fig. 3(b)] 
because the effective dipolar field is not symmetric with 
respect to the symmetry axis of the trap. 

For B ~ 50 mG, as used in the above calculations, the 
quadratic Zeeman energy (~ 0.01 nK) is much less than 
the dipolar interaction energy (edd ~ 0.2 nK), and hardly 
affects the texture formation. Using the parameters from 



the Berkeley experiment Q , we found that the results are 
qualitatively insensitive to the magnetic field for B < 
200 mG. 

The texture-formation dynamics can be observed also 
in other atomic species. We have obtained qualitatively 
the same results for 23 Na and 52 Cr atoms as those for 
87 Rb: transverse spin texture first develops due to the 
inhomogeneity of 5o, followed by a change of the local 
longitudinal magnetization. The time scale of the tex- 
ture formation depends on the strength of the dipole- 
dipole interaction; thus, for the case of 52 Cr atoms with 
a magnetic dipole moment of 6/iB, spin texture develops 
144 times faster than for spin-1 alkali atoms. 

In conclusion, we have shown that the dipole-dipole 
interaction between atomic dipole moments can induce 
observable effects on spin dynamics under an experimen- 
tally accessible magnetic field (~ 100 mG), even though 
the short-range interaction and Zeeman energy domi- 
nates the dipolar interaction. Our proposal can be veri- 
fied by taking a tilt angle of the spin other than 90° in 
the Berkeley experiment [7]. 
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